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ABSTRACT 



This thesis investigates the accuracy of a recently proposed passive bearings-only 
Target Motion Analysis (TMA) procedure. The primary method of analysis is to com- 
pare computer generated positions of a Target that is moving with a constant course and 
speed, with the procedurally derived estimated positions. 

A computer model was developed which simulated several possible interactions be- 
tween the Target and Own Ship. Estimated parameters of the Target track were com- 
puted using the procedure under analysis. These values were compared to the "true" 
values generated by the simulation. 

An analysis of the TMA procedure as originally proposed showed that it failed to 
accurately estimate the target track parameters. However with some modifications, the 
accuracy improved significantly and it is felt that the procedure can accurately estimate 
target range (but not necessarily course and speed) for some target geometries. 
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I. INTRODUCTION TO THE PASSIVE RANGING PROBLEM 



A. BRIEF HISTORY 

The importance of doing Target Motion Analysis (T.MA)i by a passive bearings- 
only tracking technique has had the attention of many tacticians over the last 45 years, 
mainly because the submarine's objective is to remain covert (undetected) during the 
tracking maneuver. This guarantees the advantage of surprise if an attack is to be per- 
formed. 

In the early part of WWII, TMA by passive bearings-only methods began to be 
used, and the following drawn from Naval Tactical Decision Aids by D. H. Wagner, 1989, 
provides a brief history of its development: 

During WWII LT. F. C. Linch developed the Linch Plot, which uses a pivotal 
relationship among bearings, bearing rate, and Target relative motion. Based on 
that, he developed a graphical method that was used until the late 1960's. It was 
abandoned because it was useful only for the short ranges contemplated by Linch 
at that time. 

During the 1950's, various human plot methods and nomographs were used for 
TMA. One of the oldest of these was the Strip Plot, later called Geographical Plot. 
This method has often been used to provide upper and lower bounds on Target 
range, based on bounds on speed. An additional much used nomographic device 
has been the Bearing Rate Slide Rule. An early TMA computer was the Position 
Keeper, which was a carryover from WWII. and used as an aid to approach and 
attack surface ships. 

In 1953, F. X. Spiess published a graphical procedure requiring four bearings and 
a maneuver by Own Ship. This method is known as Spiess Plot. It remains in use 
in surface ship TMA. 

The 1954 command thesis of LT. J. F. Fagan derived a four-bearing TMA sol- 
ution. which was a transcendental system of three equations with three unknowns. 
To reduce to computability by slide rule, he assumed that Own Ship motion during 
the first three bearings was approximately zero which was probably satisfactory for 
diesel operations. 

In 1957. LT. J. J. Ekelund devised one of the most famous TMA methods, now 
called Lkelund Ranging. 

Between 1967 and 1969. D. C. Bossard with a group of USX officers, made an 
analysis of the fundamentals of bearings-only ranging and the effect of Own Ship 
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maneuvers on ranging accuracy. Through this work and trials at sea, they observed 
that ranging errors could be eliminated by judicious choice of the time for which the 
range was estimated. Also, by selecting some particular maneuvers, this best time 
could be controlled to be in the past, present, or future; then they applied this 
method to various forms of bearings-only ranging, especially to Ekelund and Spiess 
methods. The general name for these procedures is Time Correction. 

In August 1968, a series of exercises was conducted aboard the USS Pargo to 
compare the different methods in use at that time. The Passive Ranging Manual 
evolved in three volumes from the report of the analysis of these exercises and was 
published later as an MVP. This is a basic reference today for TMA. 

In the late 1960's the CHURN method appeared, which uses a large number of 
bearing observations and regression methods to do TMA. This method was devel- 
oped by General Dynamics/ Electric Boat and Librascope. 

In 1969. H. W. Headle, J. Di Russo, and E. Messere developed the Manual 
Adaptive TMA Evaluator (MATE). This method begins with an input estimate of 
Target course, speed and range, and from there estimated bearings are computed 
and compared against observed bearings. Differences between these bearing values 
are plotted against time. If these differences are distributed roughly symmetrically 
as white noise about zero, while Own Ship changes course or speed, then the trial 
solution is a good one; otherwise the MATE operator adjusts the input Target pa- 
rameters. 

In the early I970's, the application of Kalman filtering appeared for the first time 
in TMA. It was a done by J. S. Davis of NUSC, based on the Ph. D. thesis of D. 

J. Murphy at Northeastern University. Later versions were developed by IBM. 

In the mid 1970's, a new TMA method known as FLIT was devised by ENS. L. 
Anderson which has seen considerable effective operational use. Its methodology 
remains classified. 

In the mid and later 1970's, ranging and tracking at long range on a sphere were 
developed by I). C. Bossard. J. B. Behrla. and L. K. Graves of Daniel H. Wagner, 
Associates. This was motivated by over-the-horizon targeting needs. 

In the late I970's and during the 1980's, further developments were made by using 
Kalman filtering and also stochastic differential equations. The Maneuvering Target 
Statistical Tracker (MTST) is a method devised from this basis by W. H. Barker of 
Daniel II. Wagner. Associates. 

From 19S5 to 19S7 the Generical Statistical Tracker (GST) was developed at 
Daniel II. Wagner, Associates. This is a PC version of MTST. 

B. PURPOSE OF THE THESIS 

In late 19SS. LCDR P. K. Peppe of the USS La Jolla (SSN 701) proposed a 
bearing-only Target Motion Analysis (TMA) procedure (here called the Bearings Ex- 
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trapolation Procedure) for use in submarines [Ref. 1]. It is the purpose of this thesis to 
test the accuracy of this procedure against a nonmaneuvering Target. 

C. PROBLEM DESCRIPTION 

A submarine (Own Ship) is trying to determine the range, course, and speed of a 
ship (Target) that is moving with a fixed course and speed. The data collected by Own 
Ship consists of passive bearings. Since the bearing data is obtained by a single tracker, 
a unique tracking solution may only be obtained by a maneuvering Own Ship [Ref. 2]. 
A TMA method which uses a two-leg maneuver can fulfill this requirement. 

The two-leg TMA procedure used in this study is depicted in Figure 1. Own Ship 
begins by receiving bearings to the Target each 20 seconds for 4 minutes. In Figure 1, 
/, is the time at which bearing i is received. After time t n , the bearings stop and Own 
Ship begins a course change. At time t i0 Own Ship is on the second leg and bearings 
resume. Bearings continue on the second leg, the last one being received at time / 41 . 
Estimates of Target course, speed and range will be derived from these 24 observed 
bearings. 

The central idea of the Bearings Extrapolation Procedure is that given three bearings 
observed at three times on a single leg. all future and past bearings can theoretically be 
determined. (This idea was known by Spiess [Ref. 3] and may predate him.) So, refer- 
ring to Figure 1, the estimated Target position at time t 0 can be obtained by crossing the 
true bearing (beginning at point t 0 ) with the extrapolated bearing from leg two (begin- 
ning at point r 0 '). 

The same as above can be done for all time points that have bearings between t 0 and 
/ 41 . Then regressing linearly on these estimated Target positions, the Target track pa- 
rameters can be determined. 

Although this procedure should give exact answers, in practice large errors can be 
generated. The reason for these inaccuracies is that bearings to the Target cannot be 
measured exactly, so the extrapolated bearings will also be in error. 
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Figure 1. Graphical Representation of the Bearings Extrapolation Procedure. 



Exactly how bearings that contain errors are extrapolated, with as much accuracy 
as possible, forward and backward in time is part of what this thesis examines. The most 
straightforward method (and the first to be analyzed) is simply to use linear regression 
on the ( B , , t) data, where B, is the bearing at time t. Alternate methods include re- 
gressing the (B , , t) data by using quadratic or arctangent models. Each of these proce- 
dures were tried and accuracies assessed. 
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The Spiess Graphical Method, Ekelund Ranging, and the CHURN Method are 
closely related to the procedure proposed by Peppe, and for that reason they are further 
explained in Appendix A. 

D. ASSUMPTIONS FOR THE SIMULATION 

The computer simulation used in this study uses the following assumptions: 

1. The Target moves at a constant course and speed throughout the tracking maneu- 
ver. 

2. Differences in depth for the Target and Own Ship are disregarded. 

3. Errors in position for Own Ship are disregarded. 

4. Errors in bearings to the Target are independent and identically distributed normal 
random variables with mean zero and standard deviation not larger than 1.5°. In 
addition, this standard deviation is constant for each simulation experiment. 

5. Bearings to the Target arc obtained each 20 seconds. 

6. Initial course, speed and range for the Target and initial course and speed for Own 
Ship are as described in Chapter II. 

E. TERMINOLOGY 

1. Lead-leg. An Own Ship course such that Own Ship speed across the line of sound 
is in the same direction as the Target speed across the line of sound at the start of 
the leg. 

2. Lag-leg. An Own Ship course such that Own Ship speed across the line of sound 
is in the opposite direction to the Target speed across the line of sound at the start 
of the leg. 

3. Leg length. The distance along either leg where bearings are received (does not 
include Own Ship turn). 

4. Raw bearing. An unsmoothcd bearing to the Target. 

5. Faired bearing. A bearing obtained after smoothing of the raw bearings. 

6. Extrapolated positions. Positions of Own Ship obtained by extending the first leg 
forward in time or the second leg backward in time. 

7. Extrapolated bearings. Bearings obtained by extrapolating actual bearings forward 
in time for the first leg or backward in time for the second leg. 

S. Lead-leg estimation. The estimated target track obtained by linear regression of 
estimated Target positions. These estimated Target positions are obtained from the 
intersection of faired bearings from the lead-leg and extrapolated bearings from the 
lag-leg. 

9. Lag-leg estimation. The estimated target track obtained by linear regression of 
estimated Target positions. These estimated Target positions are obtained from the 
intersection of faired bearings from the lag-leg and extrapolated bearings from the 
lead-leg. 



10. Lead angle. Angle between the line of sound and the course of Own Ship when it 
is leading the Target. 

11. Lae angle. The angle between the line of sound and the course of Own Ship when 
it is lagging the Target. 

12. Angle on the bow. The angle between the line of sound and the course of the 
Target. 
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II. TESTING THE BEARING EXTRAPOLATION PROCEDURE 



A. INITIAL CONDITIONS 

In order to perform a thorough analysis of the Bearing Extrapolation Procedure, 
5,184 separate simulation experiments were conducted. Each experiment used one of the 
possible combinations of the following parameters: 

1. The initial Target positions were all due north from Own Ship (i. e., 000°) at one 

of three different initial distances: large (= 60,000 yds), medium ( = 30,000 yds), 

and short ( = 10,000 yds). 

2. Six courses for the Target were used: 030°, 060°, 090°, 120°, 150°, and 175°. 

3. Three speeds for the Target: 25, 15, and 5 knots. 

4. Four courses for Own Ship in the lead-leg: 060°, 070°, 080°, and 090°. 

5. Four courses for Own Ship in the lag-leg: 270°, 280°, 290°, and 300°. 

6. Three speeds for Own Ship: 5, 10, and 15 knots. 

7. Two maneuvers by Own Ship: lead-lag and lag-lead. 

8. Four bearing error standard deviations: 0°, 0.5°, 1°, and 1.5°. 

It is important to mention that because of problem symmetry’ and the fact that both 
lead-lag and lag-lead maneuvers are simulated, there was no need for simulating courses 
for the Target from 1S0° to 360°. 

The randomness in the simulation was introduced only through errors in the re- 
ceived bearings. 

B. DESCRIPTION OF THE SIMULATION 

The simulation was written in Fortran 77 [Ref. 4] and executed on an IBM 3033 
mainframe. Appendix B has a complete listing of all variables of the simulation pro- 
gram. Appendix C is a listing of the simulation program code. 

For each combination of initial conditions, the simulation was repeated 100 times 
with different seeds for pseudorandom-generation of bearing error. For each of the re- 
petitions. the following steps were conducted: 

1. Target positions were simulated every 20 seconds for the entire maneuver. 

2. For the first leg the following data was generated: 

a. Own Ship position every 20 seconds. 

b. Bearing to the Target every 20 seconds (normal errors are added). 
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c. Own Ship track extended beyond the turn. 

3. For the first leg the following data was then computed: 

a. Coefficients for linear and quadratic fits for bearings versus time. 

b. Faired bearings and extrapolated bearings using the results from linear or 
quadratic fit. 

4. Own Ship course was determined for the second leg; based on the required second 
leg. lead or lag angle, and the line of sound specified by the extrapolated bearing 
at the middle of the maneuver (/ 20 ). 

5. Simulations and computations for the second leg were completed in the same 
fashion. 

6. Estimated positions of the Target in two coordinates (X and Y) were computed by 
determining the point of intersection of faired bearings with their corresponding 
extrapolated bearings. 

7. Linear and orthogonal regression is conducted on the estimated Target positions 
(in X versus Y coordinates), to determine the Target course. 

S. Range to the Target at the end of the maneuver was obtained by computing the 
distance between Own Ship's actual position at the end of the second leg and the 
estimated Target position. 

9. Finally, Target speed was obtained by computing the distance between the esti- 
mated Target position at time / 0 and the estimated Target position at the end of the 
second leg (/ 41 ), on the fitted Target track, and then dividing this distance by 
hi h' 

The simulation of the bearing errors was done using the Linear Congruential 
Method for uniform pseudorandom numbers, combined with the Box and Muller 
Method to produce normal variates (Refs. 5,6]. These procedures were selected to allow 
the simulation to be performed on a personal computer. 

C. TESTING THE BEARING EXTRAPOLATION PROCEDURE UNDER IDEAL 
CONDITIONS 

Here bearings to the Target were considered without error. Also the bearing rates 
immediately before the turn (used to extrapolate bearings forward in time) and imme- 
diately after the turn (used to extrapolate bearings backward in time) were computed 
exactly with no errors introduced. The only source of error in this case was the linear 
bearing extrapolation. This case is intended to be an optimistic bound on the accuracy 
of the bearing extrapolation procedure when linear regression is used. 

The measure of effectiveness (MOE) for the accuracy of this procedure was the 
percentage of times that the estimated error was within 10% for range. 20° for course, 
and 10? o for speed from the actual parameters of the Target track. These results are 
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presented in Table 1 and Table 2 for lead-lag and lag-lead maneuvers respectively. 
Each number in these Tables resulted from the examination of 18 simulation exper- 
iments (six Target courses and three Target speeds). For example, when the initial dis- 
tance was large, lead and lag angles 70°, and Own Ship speed 15 knts, then five of the 
18 simulations experiments (or 27.8%) resulted in the calculated Target speed being 
within 10% of the actual Target speed. Note that without bearing error, the simulation 
experiments are deterministic. 
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Table 1. ACCURACY IN IDEAL CONDITIONS FOR BEARINGS EXTRAPO- 
LATION PROCEDURE (LEAD-LAG MANEUV ER) 



Own Ship 
speed 
(knts) 


Lead and 
lag angle 
(°) 


Initial dis- 
tance 


Percentage of times the absolute error is 
smaller than 10% of actual for range and 
speed or 20° of actual for course. 


Ranee 

(10%) 


Course 

(20°) 


Speed 

(10%) 






Large 


100 


50.0 


27.8 




70 


Medium 


100 


50.0 


27.8 






Short 


38.9 


50.0 


22.2 






Large 


100 


50.0 


33.3 


15 


SO 


Medium 


100 


55.6 


33.3 






Short 


50.0 


50.0 


27. S 






Large 


100 


55.6 


27.8 




90 


Medium 


100 


55.6 


27.8 






Short 


50.0 


66.7 


33.3 






Large 


100 


44.4 


44.4 




70 


Medium 


94.4 


38. 9 


33.3 






Short 


55.6 


38.9 


22.2 






Large 


100 


44.4 


33.3 


10 


80 


Medium 


94.4 


44.4 


33.3 






Short 


55.6 


44.4 


27.8 






Large 


100 


55.6 


j j.3 




90 


Medium 


94.4 


50.0 








Short 


61.1 


50.0 


27. S 






Large 


S3. 3 


44.4 


22.2 




70 


Medium 


72.2 


44.4 


27.8 






Short 


44.4 


44.4 


27. S 






Large 


SS.9 


33.3 


27.8 


5 


SO 


Medium 


77. S 


33.3 


27.8 






Short 


44.4 


44.4 


2 ~> .2 






Large 


SS.9 


33.3 


27. S 




90 


Medium 


77.8 


33.3 


22.2 






Short 


44.4 


33.3 


16.7 
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Table 2. ACCURACY IN IDEAL CONDITIONS FOR BEARINGS EXTRAPO- 
LATION PROCEDURE (LAG-LEAD MANEUVER) 



Own Ship 
speed 
(knts) 


Lead and 
lac ancle 

(°f 


Initial dis- 
tance 


Percentage of times the absolute error is 
smaller than 10% of actual for range and 
speed or 20° of actual for course. 


Ranee 

(10%) 


Course 

(20°) 


Speed 

(10%) 






Large 


SS.9 


77.8 


27. S 




70 


Medium 


66.7 


77.8 


27.8 






Short 


50.0 


50.0 


33.3 






Large 


SS.9 


94.4 


3S.9 


15 


SO 


Medium 


72.2 


SS.9 


50.0 






Short 


50.0 


33.3 


33.3 






Large 


94.4 


94.4 


44.4 




90 


Medium 


72.2 


S8.9 


44.4 






Short 


3S.9 


21 .% 


21 .% 






Large 


83.3 


88.9 


22.2 




70 


Medium 


66.7 


94.4 


22.2 






Short 


44.4 


50.0 


21 .% 






Large 


S3. 3 


94.4 


jj. 3 


10 


SO 


Medium 


61.1 


SS.S 








Short 


50.0 


44.4 


22.2 






Large 


77. S 


94.4 


3S.9 




90 


Medium 


61.1 


SS.9 


3S.9 






Short 


50.0 


3S.9 


16.7 






Large 


72.2 


11 .% 


33.3 




70 


Medium 


50.0 


11.% 


27.8 






Short 


3S.9 


66.7 


27. S 






Large 


72.2 


11 .% 


33.3 


5 


SO 


Medium 


50.0 


11 .% 


21 .% 






Short 


44.4 


66.7 


27. S 






Large 


72.2 


77.8 


3S.9 




90 


Medium 


50.0 


11 .% 








Short 


44.4 


66.7 


22.2 



11 



As can be seen in Table 1 and Table 2, the following results were obtained under 
ideal conditions: 

1. Range estimation is good, but only for medium or large initial distances. 

2. Accuracy improves with increases in Own Ship speed. 

3. Course and speed estimation is less accurate than range estimation. 

4. Accuracy improves with increases in lead and lag angle. 

5. Lead-lag maneuvers give better range estimates than do lag-lead, but not neces- 
sarily better course and speed estimates. 

It is felt that the procedure performed poorly for short initial distances, in part be- 
cause the bearing versus time curve is more nonlinear in this case. 

D. TESTING WITH BEARING ERRORS 

For these experiments, the program was run with bearing errors added to the simu- 
lated bearings. Also the bearing rate for each leg was computed using a linear least 
squares technique. 

Only lead or lag angles of 90° were considered, since it was determined in the pre- 
vious step that this is the best course for Own Ship to follow. Accuracy in course and 
speed were not considered since the errors for these parameters were too large, even in 
ideal conditions. An Own Ship speed of 5 knots was not considered for the same reason. 

The MOE for these experiments was the percentage of times (runs) that the esti- 
mated range at time / 4 , was within 5%, 10%. and 20% of the actual range. Each number 
in Table 3 and Table 4 results from 1,800 trials (six Target courses, three Target speeds 
and 100 replications per (course, speed) combination). 
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Table 3. ACCURACY IN REAL CONDITIONS FOR BEARING EXTRAPO- 
LATION PROCEDURE (LEAD-LAG MANEUVER) 



Own Ship 
speed 
(knts) 


St. Dev. 
bearing er- 
ror 


Initial dis- 
tance 


Percentage of times the absolute error is 
smaller than 5%, 10%, and 20% of ac- 
tual for range. 


5% 


10% 


20% 






large 


6.2 


11.8 


25.5 




1.5° 


medium 


12.3 


24.4 


45.4 






short 


15.9 


30.6 


53.7 






large 


S.9 


18.4 


36.6 


15 


1.0° 


medium 


15. S 


37.8 


60.3 






short 


18.7 


38.3 


60.4 






large 


16.9 


33.7 


61.9 




0.5° 


medium 


26.5 


52.5 


82.2 






short 


19.4 


39.2 


62.8 






large 


4.2 


8.8 


18.2 




1.5° 


medium 


8.4 


17.0 


33.9 






short 


12.6 


24.4 


51.2 






large 


5.8 


12.3 


26.4 


10 


1.0° 


medium 


10.7 


23.5 


45.3 






short 


15.1 


30.2 


59.7 






large 


10.8 


22.9 


45. S 




0.5° 


medium 


19.0 


38.9 


68.7 






short 


20.3 


38.8 


69.4 
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Table 4. ACCURACY IN REAL CONDITIONS FOR BEARING EXTRAPO- 
LATION PROCEDURE (LAG-LEAD MANEUVER) 



Own Ship 
speed 
(knts) 


St. Dev. 
bearing er- 
ror 


Initial dis- 
tance 


Percentage of times the absolute error is 
smaller than 5%. 10%, and 20% of ac- 
tual for range. 


5 °/ 

~Vo 


10% 


20% 






large 


4.7 


1 1.3 


23.2 




1.5° 


medium 


10.3 


19.7 


36.5 






short 


6.7 


12.3 


23.9 






large 


7.11 


15.9 


34.3 


15 


1.0° 


medium 


11.7 


24.2 


46.6 






short 


6.7 


12.2 


22.2 






large 


16.2 


31.8 


56.1 




0.5° 


medium 


17.1 


30.9 


54.9 






short 


5.3 


11.6 


21.4 






large 


4.0 


7.4 


18.1 




1.5° 


medium 


7.3 


13.7 


26.9 






short 


5.7 


11.7 


21.9 






large 


5.4 


10.8 


23.4 


10 


1.0° 


medium 


9.0 


17.3 


35.9 






short 


5.0 


10.S 


23.2 






large 


11.0 


22.2 


42.9 




0.5° 


medium 


12.4 


25.1 


46.2 






short 


3.S 


9.7 


23.6 



From Table 3 and Table 4, the following results are obtained: 

1. The higher the Own Ship speed, the better the accuracy. 

2. The best accuracy is obtained for medium initial distances. 

3. The smaller the bearing error, the better the accuracy. 

4. A lead-lag maneuver gives better accuracy than lag-lead maneuver. 

It is noted that orthogonal regression was also tried on the estimated Target posi- 
tions, but with results similar to those obtained with linear regression. Subsequent ex- 
periments used only linear regression. 
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E. ANALYSIS OF THE RESULTS 



It appears that even in the absence of bearing errors, the Bearing Extrapolation 
Procedure gives poor estimates of Target course and speed. And when bearing errors 
are introduced, the range errors also become unacceptable. It is the purpose of this 
section to determine why these poor results were obtained. 

1. Differences Between Lead-Leg and Lag-Leg-Estimation 

Figures 2 to 5 show actual and estimated Target positions for four encounter 
geometries with no error in bearing. So the only errors introduced are due to the linear 
extrapolation of bearings. It is clear that lag-leg estimation is much more accurate than 
the lead-leg-estimation. It is also clear that any attempt to regress estimated Target 
positions over both legs can result in large errors. 
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Figure 2. Actual (.) and Estimated (*) Target Track for Lead- Lag Maneuver. 
(Target Course = 060°) 
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Figure 3. Actual (.) and Estimated (*) Target Track for Lead-Lag Maneuver. 
(Target Course = 150°) 
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Figure 4. Actual (.) and Estimated (*) Target Track for Lag-Lead Maneuver. 
(Target Course = 030°) 
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Figure 5. Actual (.) and Estimated (*) Target Track for Lag-Lead Maneuver. 
(Target Course = 120°) 



2. Explanation of Why Lag-Leg Estimation is Better than Lead-Leg Estimation 

To help explain why the lag-leg estimation was better than lead-leg estimation, 
scatter plots of bearings versus time for numerous encounter geometries were obtained. 
Figures 6, 7, and 8 are typical. In Figure 6, a lead-lag encounter is represented. It was 
observed that during the lead-leg (first leg), the bearing rate was generally smaller than 
that of the lag-leg (second leg). Also, and perhaps more importantly, the change in 
bearing rate (bearing acceleration) was also smaller in the lead-leg. When this happens, 
it is possible to use linear regression to extrapolate the lead-leg bearings accurately into 
the future. This makes the lag-leg estimation more accurate. It is clear from Figure 6 
that if linear regression were used to extrapolate the lag-leg bearings into the past, sig- 
nificant errors would result, and these errors would likely cause the lead-leg estimation 
to be less accurate. 

In the same fashion, Figure 7 shows a plot of bearings versus time for a typical 
lag-lead encounter. Here also, the lead-leg (second leg) has a smaller bearing rate and 
bearing acceleration. So linear regression can be used to extrapolate accurately the 
lead-leg bearings backward in time, again making the lag-leg estimation most accurate. 



Unfortunately, however, there are encounter geometries where neither the lead 
nor the lag-legs show small bearing acceleration. For example, Figure 8 shows a case 
where the Target is at high speed, short initial distance, and on course 160°. Here both 
lead and lag-leg bearings are nonlinear, making the Bearing Extrapolation Procedure 
unusable. Extensive experimentation indicates that when all the following conditions 
hold: 

Target course > 150 °, 

Target speed > 5 knots, 

Target range > 40,000 yds; 

then the lead-leg bearings (and the lag-leg bearings) are not significantly linear to be 
accurately extrapolated, and as a result, lag-leg estimation (as well as lead-leg esti- 
mation) is poor. 

Finally, it is noted that for very long initial distances (> 40,000 yds), both the 
lead and lag-legs have approximately a linear change in the bearings, making both lead 
and lag-leg estimations reasonably accurate. 




Figure 6. Bearings in Both Legs for Lead-Lag Maneuver (Target course = 060°) 
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Figure 7. Bearings in Both Legs for Lag-Lead Maneuver (Target Course = 140°) 




Figure 8. Bearings in Both Legs for Lead-lag Maneuver at Short Distance and 
Closing Target. 
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F. CONCLUSIONS 



1. Even when bearings were without error, the Bearing Extrapolation Procedure was 
unsuccessful in estimating Target course and speed. However, range estimates were 
reasonable. 

2. When bearing errors were introduced, range estimates also became inaccurate. 

3. For most encounter geometries, lag-leg estimation gave more accurate results than 
did lead-leg estimation. However, for some important cases both lag and lead-leg 
estimations performed poorly. 
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III. POSSIBLE IMPROVEMENTS, AND ANALYSIS OF THEIR 

ACCURACIES 



The lack of accuracy of the proposed procedure, even under ideal conditions, led to 
attempts to improve its performance. The following are some suggested minor changes 
and two suggested major changes that were found to improve the accuracy of the Bear- 
ing Extrapolation Procedure. 

A. MINOR CHANGES 

1. Increase in Leg Length 

In an attempt to improve performance, several different leg lengths were tried. 
It was discovered that leg lengths of 8 to 12 minutes gave the best results. Shorter legs 
gave too little data and longer legs resulted in nonlinear changes in lead-leg bearings. 
Based on these tests, ten minutes was selected as the leg length for subsequent simu- 
lation experiments. 

2. Use of a Quadratic Model to Fair the Lag-Leg 

Another important change from the original simulation that was tried was to 
use a quadratic model to fair the bearings in the lag-leg. This produced better accuracy 
in the estimated Target track parameters. It is noted that current manual TMA methods 
(in particular, the Time-Bearing Plot) use linear fairing. 

3. Increase in Turn Time 

Contrary to Ekelund ranging, if an instantaneous turn is considered, large errors 
in the estimated Target track parameters are generated. This is because some actual and 
corresponding extrapolated bearings are too close to each other (i. e., a small baseline). 
The intersection of those bearings often results in a estimated Target position far away 
from the true position of the Target. 

Based on these considerations, it was determined that the time required for Own 
Ship to turn should be increased. After many trials, it was found that the best turn du- 
ration is ten minutes, if leg length is also ten minutes. For longer times, the bearing rate 
in the lead-leg does not remain low and constant, reducing the accuracy of the lead-leg 
bearing extrapolation. 



B. ELIMINATION OF THE LEAD-LEG-ESTIMATION 

As shown in Chapter II, when linear extrapolation of the bearings from the lag-leg 
is used, increased errors are generated in the estimation of the Target track parameters. 
This occurs because the lag-leg generates a high and variable bearing rate compared to 
the lead-leg. A simple solution to the problem is to drop the estimated Target positions 
generated by the extrapolation of the lag-leg bearings. 

One problem with this idea is that a lead-lag maneuver will then produce more ac- 
curate estimates of Target position at the end of the maneuver than will a lag-lead ma- 
neuver. This occurs since in the lead-lag case, the estimated Target positions are 
obtained for times during the second (i.e., the lag) leg. And in the lag-lead case, the es- 
timated Target positions are for times in the first (again, the lag) leg. So to use a lag-lead 
maneuver to estimate Target positions at the end of the maneuver requires that these 
positions be extrapolated in time using derived estimates of Target course and speed. 
This leads to accumulated errors in the final Target position. 

To test for the possible improvements that might result when using only lag-leg es- 
timation. simulation experiments were conducted with the following parameters: 

1. Speed of Own Ship: 5, 10, and 15 knots. 

2. Speed cf the Target: 5, 15, and 25 knots. 

3. Initial range: large ( = 60.000), medium ( = 30,000), and short ( = 10,000) yds. 

4. Only the lead-lag maneuver is considered. 

5. Only the initial lead angle of 90° is cosidered. 

6. Leg length: 4 minutes. 

7. Time between legs: 6 minutes. 

S. Standard deviation bearing error: 0°, 0.5°. 1.0°, and 1.5°. 

9. Bearings for the first leg were faired and extrapolated linearly. 

10. Bearings in the second leg were faired linearly. 

It is important to note that other values for lead angle were not considered because, 
as discussed in Chapter II, it was found that the best results were obtained for a 90° lead 
angle. 

In order to show how the minor changes improve the performance of this major 
suggested change (i. e., elimination of the lead-leg estimation), the simulation was run 
for three dilferent cases: 



1. Without Minor Changes 

The simulation program was run 100 times for each combination of the above 
parameter values. It was determined that for the ideal case (0° bearing errors) the ac- 
curacy of the predicted Target parameters was adequate, but with bearing errors intro- 
duced, large errors in all three Target parameters were produced. So it was concluded 
that simply eliminating the lead-leg estimation was not a sufficient improvement. 

2. With Increased Leg Length and Lag- Leg Faired with a Quadratic Model 

Here leg length was increased to ten minutes and the lag-leg was faired with a 
quadratic model. The same simulation experiments were run as before. Only the results 
for lead-lag maneuvers are presented in Table 5. Lag-lead maneuvers gave uniformly 
worse results for reasons discussed above. Additionally, and more surprisingly, the ac- 
curacy for the lag-lead maneuver was not as good as the lead-lag, even when estimating 
the Target parameters at the end of the first leg. This was a confirmation that, in gen- 
eral, the lead-lag maneuver gives better results than do lag-lead [Ref. 7: pp. 3]. 

Results for estimated Target course and speed were very precise only when the 
bearing error was zero and when Own Ship speed was greater than ten knots. When 
errors in bearings were included, the accuracy in estimated Target course and speed was 
tremendously reduced. For example, for large initial distance and 15 knots Own Ship 
speed the following was obtained: 

1. With no bearing error, the error in estimated speed was always less than 10% and 
the error in course was less than 10° in 83% of the different geometries. 

2. With a small bearing error (standard deviation = 0.5°), only 33.4% of the runs had 
an error in speed of less than 10% . In 15.1% of the runs, the error in course was 
less than 10°. 

One of the conclusions of this study is that the Bearing Extrapolation Procedure, even 
as modified here, does not produce accurate estimates of Target course and speed. With 
Target range, however, more success was obtained. Table 5 shows the percentage of 
simulation runs that produced an estimated Target range at the end of the maneuver 
within 5'%. 10%, and 20% of the true range. 
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Table 5. ACCURACY IN ESTIMATED RANGE WHEN ONLY 
LAG-LEG-ESTIMATION IS USED. (LEG LENGTH - 10 MIN). 



Own Ship 
speed 
(knts) 


St. Dev. 
bearing er- 
ror 


Initial dis- 
tance 


Percentage of times the absolute error is 
smaller than 5%, 10%, and 20% of ac- 
tual for range. 




5% 


10% 


20% 






large 


21.6 


41.7 


72.1 




1.5° 


medium 


2S.3 


54.6 


82.2 






short 


10.9 


19.7 


40.9 






large 


30.6 


57.1 


84.6 




1.0° 


medium 


35.4 


62.9 


86.6 


15 




short 


10.7 


18.6 


40.2 




large 


49.9 


78.8 


93.7 




0.5° 


medium 


43. S 


71.1 


88.6 






short 


10.3 


IS. 9 


39.6 






large 


72.2 


88.9 


94.4 




0.0° 


medium 


55.6 


77.8 


8S.9 






short 


n.i 


22.2 


38.9 






large 


15.2 


30.6 


55.9 




1.5° 


medium 


20.2 


39.4 


70.7 






short 


16.4 


32.4 


5 1 .3 






large 


21.0 


41.2 


70.1 




1.0° 


medium 


26.2 


49.5 


76.6 


10 




short 


16.1 


32.7 


52.6 




large 


35.4 


62.1 


88.7 




0.5° 


medium 


35.4 


61.3 


79.5 






short 


15.3 


32.3 


52.5 






large 


72.2 


77. S 


94.4 




0.0° 


medium 


3S.9 


72.2 


77.8 






short 


16.7 


27.S 


50.0 



From Table 5 it can be seen that the larger the Own Ship speed , the better is 
the estimate of range. This is shown in Figure 9 which is the plot of the accuracy in 
estimated range versus standard deviation of bearing errors when the Own Ship speed 
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is 10 and 15 knots. This accuracy is the percentage of times that the absolute error is 
within 20% of the actual for range at the end of the maneuver. 




Figure 9. Comparison of Accuracy for Different Own Ship Speeds by Using only 
the Lag-Leg-Estimation. 



Figure 10 shows the comparison of the results of these suggested improvements with the 
results obtained in Chapter II for the central idea. It appears that for short initial dis- 
tances the accuracy of the procedure decreases, but for large and medium distances, even 
with large bearing error (standard deviation = 1.5°), the accuracy increases consider- 
ably. Note that the accuracy for short distances was not adequate in either case. 
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Figure 10. Comparison of the Accuracy Between Results by Using only the Lag- 
Leg-Estimation and the Results by Using the Central Idea. 



3. With Previous Minor Changes Plus an Increase in Turn Duration 

In addition to the two minor changes made before, here the turn duration was 
increased to ten minutes and the simulation experiments were repeated. Table 6 shows 
the results for this case. 
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Table 6. ACCURACY IN ESTIMATED RANGE WHEN ONLY 
LAG-LEG-ESTIMATION IS USED. (LEG LENGTH = 10 MIN. 
TURN TIME = 10 MIN.) 



Own Ship 
speed 
(knts) 


St. Dev. 
bearing er- 
ror 


Initial dis- 
tance 


Percentage of times the absolute error is 
smaller than 5%, 10%, and 20% of ac- 
tual for range. 


5% 


10% 


20% 






large 


25.5 


47.7 


7S.3 




1.5° 


medium 


30.6 


54.3 


83.3 






short 


12.1 


21.7 


38.6 






large 


34.8 


62.0 


87.8 




1.0° 


medium 


37.3 


59.9 


86.7 


15 




short 


13.3 


22.9 


38.2 




large 


44.3 


73.5 


91.3 




0.5° 


medium 


40.7 


62.1 


87.8 






short 


14.4 


24.2 


3S.6 






large 


66.7 


88.9 


94.4 




0.0° 


medium 


50.0 


61.1 


8S.9 






short 


16.7 


22.2 


38.9 






large 


1 S.3 


36.2 


65.1 




1.5° 


medium 


24.6 


46.2 


72.0 






short 


14.7 


2S.7 


50.3 






large 


25.9 


4S.2 


7S.4 




1.0° 


medium 


29.6 


53.9 


75.5 


10 




short 


13.3 


28.1 


50.6 




large 


31. S 


59.0 


85.5 




0.5° 


medium 


32.2 


5S.6 


76.9 






short 


13.2 


26.6 


50.9 






large 


66.7 


77.8 


94.4 




0.0° 


medium 


s -> 


66.7 


77. S 






short 


16.7 


22.2 


44.4 



By comparing the result's obtained for the accuracy in range in Table 5 (Own 
ship turn = 6 minutes) and Table 6 (Own ship turn = 10 minutes), it is appears that 
small improvements in precision are obtained by using ten minutes for the turn. 
Figure 1 1 is a graphic representation of this fact. 



OWN SUP SPEED - IS Knt*. TURN TIME - • Un» 
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Figure 11. Comparison of Accuracy for Different Own Ship Turn Times by Using 
only the Lag-Leg-Estimation. 

a. Effect of Target Courses Greater than 15(T when Lag-Leg Estimation and 
Minor Changes are Used 

It was stated in Chapter II that the bearing rate for the lead-legs does not 
remain constant when the Target course is greater than 150° (initial angle on the bow 
less than 30°). This is because range during the maneuver decreases to the point where 
the lead-leg bearings do not maintain linearity. To determine the effect on the esti- 
mation of Target range that this behavior produces, simulations were conducted as be- 
fore except that results were obtained for 12 separate Target courses between 030° and 
175°. Previous results were averaged over Target course, so the effect of a particular 
Target course was not observed. 

The results of one simulation experiment are depicted in Figure 12 and 
presented in Table 7. The parameters used to obtain these results were the following: 

1. Initial distance: 30,000 yds. 

2. Target couses: 030°, 050°, 070°, 090°, 110°, 130°, 150°, 155°, 160°, 165°, 170°, and 
175°. 

3. Target speeds: 25, 15, and 5 knots. 
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4. Own Ship speed: 15 knots. 

5. Standard deviation of bearing error: 1.5°. 

6. Leg length: 10 minutes. 

7. Turn time: 10 minutes. 




Figure 12. Fraction of Times Absolute Error is within 20% of Actual for Range, for 
Different Target Courses and Speeds. 
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Table 7. FRACTION OF TIMES ABSOLUTE ERROR IS WITHIN 20% OF 
ACTUAL FOR RANGE, FOR DIFFERENT TARGET COURSES AND 
SPEEDS. 



Target 

speed 

(knts) 


Target course (°) 


030 


050 


070 


090 


110 


130 


150 


155 


160 


165 


170 


175 


25 


0.87 


0.83 


0.72 


0.98 


1.00 


1.00 


0.78 


0.62 


0.00 


0.00 


0.00 


0.00 


15 


0.86 


0.95 


0.95 


0.97 


0.99 


1.00 


0.94 


0.62 


0.38 


0.26 


0.13 


0.05 


5 


0.94 


0.96 


0.94 


0.97 


0.95 


0.93 


0.93 


0.94 


0.91 


0.97 


0.91 


0.95 



From Figure 12 and Table 7. it can be seen that for a Target speed of 5 
knots, approximately equal accuracy was obtained for all Target courses. But for Target 
speeds of 15 and 25 knots, much less accuracy was obtained when the Target course was 
between 150° and 175°. For Target courses less than 150°, there was less observed 
change in the estimation of Target range. 

C. USE OF AN ARCTANGENT MODEL TO EXTRAPOLATE THE LAG-LEG 
BEARINGS 

As presented so far. a major problem with the Bearing Extrapolation procedure is 
that the lag-leg bearings cannot be linearly extrapolated forward or backward in time 
while maintaining accuracy. The solution to this problem suggested in the previous se- 
veral sections was simply not to use extrapolated bearings from the lag-leg. The problem 
with this approach is that tracking information is lost when these extrapolated bearings 
are not used. In this section, an attempt will be made to recover these extrapolated 
bearings by using an arctangent, rather than linear, extrapolation model. 

It is known that without change in Own Ship course, bearing to a constant course 
and speed moving Target is of the form 



B t = Dq + tan [ 



R 






cpa 



(3.0) 



where 

B 0 is the bearing to the Target at CPA2, 
- Closest Point of Approach 
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Y 0 is the relative speed between Target and Own Ship, 

R Cfa is the range to the Target at CPA, 

t Q is the time at which CPA is to occur. [Ref. 8: pp. 2] 



Assuming that the extrapolated bearings for the lag-leg follow this arctangent 
model, two different attempts were made in order to include these bearings in the sol- 
ution of the problem. 

I. Fitting the Arctangent Model by Using Least Squares 

An attempt was made to use a least squares technique to find the values of B 0 , 



this approach was ultimately unsuccessful, there may be another least squares procedure 
that will prove more robust. To help any subsequent researcher, the details of this at- 
tempt are presented to illustrate what has already been examined. 



R, 



and r 0 in (3.0) which produced the best fit to the observed bearing data. Although 



‘cpc 



Let 



o 



R. 



xpa 



Now make the change of variable 



T,= tan(Z?,- B 0 ), 



or 



y, = a(r - t 0 ). 



Then applying least squares, it is necssary to minimize 



n 



ss = 




' ( — ca + at 0 ) 2 . 



A necessary condition is 



n n 



= -:Y( V, - at + a t 0 )t + 2 Y( Y, - ai + at 0 )i Q = 0, (3.1) 

CO. 1 — J 1 — J 




and 
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(3.2) 



dSS 
oi o 



n 

= 2 'Yji Y, - at + at 0 )a = 0. 

f=i 



Since a cannot be zero unless relative speed (F 0 ) is zero, from (3.2), 



n 

^{Y, - at + at 0 ) = 0. 

?=i 



Dividing (3.3) by n and solving for a, 



Y 




(3.3) 



(3.4) 



Now substituting (3.3) in (3.1) gives 




(3.5) 



Substituting (3.4) in (3.5) and solving for t 0 yields 



n 




^0 n n 

T ,y - Y T,' 

1 = 1 1 



(3.6) 



Because the determination of Y, for all /. requires knowledge of the value of the 
bearing at CPA (Z? 0 ). an iterative procedure was attempted where the first value for B 0 
was guessed, then with equations (3.5) and (3.6), a and r 0 were computed. Then bv using 



n 

= 7T - tan -1 0(/ - r 0 )]}, 

;=i 
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a new value for B 0 was obtained. This iteration was continued in an attempt to converge 
on the true B 0 . Unfortunately for even small errors in bearings, (standard deviation = 
0.5°), the procedure would not reliably converge. 

2. Estimation of the Parameters for the Arctangent Model by Solving a System of 
Linear Equations 

It is known that with three different bearings to a Target that is moving with a 
constant course and speed, the unknown parameters of the equation 

B, ~ Bo + tan ' C a (t — to)] 

can be determined. [Ref. 7: pp. 2] 

Now let B u B 2 , and Z? 3 be the bearings obtained at times r,, i 2 , and i } respec- 
tively. Then it can be written 

= B 0 + tan -1 [a(r, - r 0 )], 



B 2 = B 0 + tan '[>(/ 2 - r 0 )], 



^3 “ ^0 + tail [ a (u — t 0 )]. 

An attempt will be made to solve for B 0 , a, and t 0 . Substracting. 



B 2 - B\= tan '[>(/ 2 - r 0 )] - tan ’[«(/, - r 0 )], 



B 3 ~ B\= tan ’[afo - r 0 )] - tan '[>.(/, - r 0 )], 



B 3 - B 2 = tan ’[a(/ 3 - / 0 )] - tan 1 [a(/ 2 - r 0 )]. 



Now taking the tangent of both sides and knowing the trigonometric identity for the 
tangent of the difference of two angles, the following can be obtained: 



tan(Z? 2 — Z?j) = 



«(*2 ~ b ) 

1 + a (U — to)(h — i 0 ) 



tan(Z? 3 — B ,) = 



a(h - b) 

1 + v 2 (h ~ l o)( r \ ~ t 0 ) 



tan(# 3 - B 2 ) = : . 

1 + a (/ 3 — to)(t 2 — to) 
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